Charge and Spin Currents of the ID Hubbard Model at Finite 

Energy 



N. M. R. Peres 1 , P. D. Sacramento 2 , and J. M. P. Carmelo 1 
Department of Physics, University of Evora, Apartado 94, P-7001 Evora Codex, Portugal 
2 Departamento de Fisica and CFIF, Instituto Superior Tecnico, Av. Rovisco Pais, P-1096 Lisboa 

Codex, Portugal 
(September 1997) 

Abstract 

The transport of charge and spin at finite energies is studied for the Hubbard 
chain in a magnetic field by means of the pseudoparticle perturbation theory. 
In the general case, this involves the solution of an infinite set of Bethe-ansatz 
equations with a flux. Our results refer to all densities and magnetizations. 
We express the charge and the spin-diffusion currents in terms of elementary 
currents associated with the charge and spin carriers. We show that these are 
the a, 7 pseudoparticles (with a = c,s and 7 = 0,1, 2, 3, ...) and we find their 
couplings to charge and spin. We also study the ratios of the pseudoparticle 
charge and spin transport masses over the corresponding static mass. These 
ratios provide valuable information on the effects of electronic correlations in 
the transport properties of the quantum system. We show that the transport 
of charge and spin in the Hubbard chain can, alternatively, be described by 
means of pseudoparticle kinetic equations. This follows from the occurrence 
of only forward-scattering pseudoparticle interactions at all energies. 

PACS numbers:71.10. Pm, 05.30. Fk,72.90.+y,03.65. Ca 
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I. INTRODUCTION 



The transport properties of strongly correlated electron models for low-dimensional con- 
ductors has been a subject of experimental and theoretical interest for over twenty years. 
Low-dimensional conductors show large deviations in their transport properties from the 
usual single-particle description. This suggests that electronic correlations might play an 
important role in these systems M, even if they are small M. Solvable one-dimensional 
many-electron models such as the Hubbard chain are often used as an approximation for the 
study of the properties of quasi-one- dimensional conductors P],H]. Although the Hubbard 
chain has been diagonalized long ago HQ, the involved form of the Bethe-ansatz (BA) wave 
function has prevented the calculation of dynamic response functions, these including the 
charge-charge and spin-spin response functions and their associate conductivity spectra. 

Information on low-energy expressions for correlation functions can be obtained by com- 
bining BA with conformal-field theory ||. On the other hand, several approaches using 
perturbation theory |J, bosonization |7||§, the pseudoparticle formalism ||, scaling meth- 
ods []10|1 , and spin- wave theory (Tl[| have been used to investigate the low-energy transport 
properties of the model away from half filling and at the metal - insulator transition [[3]. 
Unfortunately, only limited information on the transport properties at finite energies has 
been obtained by numerical methods [|T2T]il[j. 



Recently, a pseudoparticle description of all the BA Hamiltonian eigenstates |L5[ has 
allowed the evaluation of analytical expressions for correlation functions at finite energy 
From these results one can obtain expressions for the absorption band edges of the 



frequency-dependent electronic conductivity, o~(lu) fll7 |. The pseudoparticle theory of Ref. 
|T5| introduces new branches of pseudoparticles and generalizes for all energy scales pre- 



vious low-energy studies JT8[. The new pseudoparticle branches are associated with heavy 
pseudoparticles. These are the quantum objects needed for the description of Hamiltonian 
eigenstates showing an energy gap relatively to the ground state. (This justifies why they 
are called heavy.) As in the case of the low-energy properties of the Hubbard model [|IBf , 
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it is possible to write the Hamiltonian in terms of a set of anticommuting pseudoparticle 
operators. Importantly, all the model eigenstates can be generated from the SO (4) ground 
state |L]|. 

In this paper we use the above generalized pseudoparticle theory to study the charge 
and spin currents of the Hubbard chain at finite energy. For that we solve the BA equations 
with a twist angle for all densities and magnetizations. Moreover, we express the charge 
and the spin-diffusion currents in terms of the elementary currents of the charge and spin 
carriers. It is shown that the latter carriers are the a, 7 pseudoparticles of the pseudoparticle- 
perturbation theory (PPT) fl5f . We evaluate their couplings to charge and spin and define 
the charge and spin pseudoparticle transport masses. The ratios of these masses over the cor- 
responding static mass provide important information on the role of electronic correlations 
in the transport of charge and spin in the ID quantum liquid. Furthermore, we find that the 
transport of charge and spin can be described by means of pseudoparticle kinetic equations. 
Our results are a generalization to finite energies of the low-energy results on transport of 
charge and spin presented in Ref . || . This is possible by means of the generalized pseudopar- 
ticle representation introduced in Ref. |T5| which is a generalization to finite-energy scales 
of the usual low-energy operator representation [H| in terms of pseudoparticles PD|-[21" . 



The paper is organized as follows. In Sec. || we present the PPT description of the 



Hilbert space associated with the zero-temperature transport of charge and spin. In Sec. [HI 
we introduce the general BA equations with a spin dependent twist angle, (f> a , valid for all 
energy scales, and solve these equations by use of the PPT. The pseudoparticle transport 
and static masses and the transport of charge and spin are studied and discussed in Sec. [IV|. 
In Sec. [V| we show that the transport of charge and spin can be described by means of kinetic 
equations. Finally, the concluding remarks are presented in Sec. [VT[ 
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II. PSEUDOPARTICLE PERTURBATION THEORY 



Let us consider the Hubbard-chain Hamiltonian with iV electrons in a magnetic field H 
and with chemical potential /i 



H = -^EtsWi- + h A + ^Et%T - 1/2] [n i4 - 1/2] 



,t 



where ct (cj CT ) creates (annihilates) one electron with spin a (here and when used as operator 
index, a =f , |, and a = ±1 otherwise), n 3jCT = Cj a Cj a is the number operator at site j, N a is 
the number of sites of the chain (since we are using periodic boundary conditions it is rather 
a ring), and CN a+ i a = C\ a . In general, we use units such that h = 1, the lattice spacing a = 1, 
and the electron charge — e = 1. The form of the interaction term accounts for the particle 
- hole symmetry of the model at half filling P,^5|. For simplicity, we consider electronic 
and magnetization densities in the domains < n < 1 and < m < n, respectively, where 
n = N/N a = [iV T + N^/Na, m = [JV T - N^/Na, and N a is the number of a electrons. 



Although the pseudoparticle description refers to all Hamiltonian eigenstates |15j, in this 
paper we restrict our study to the Hilbert subspace involved in the zero-temperature charge 
and spin frequency-dependent conductivity ||17|| . This is spanned by all the Hamiltonian 
eigenstates contained in the states j^\GS), where \GS) denotes the ground state and are 
the charge (( = p) and spin (£ = a z ) current operators (given by Eqs. (|27| ) and ( pE|) below, 
respectively). Since these current operators commute with the six generators of the 77-spin 
and spin algebras [^(| , our Hilbert subspace is in the present parameter space spanned only 
by the lowest-weight states (LWS's) of these algebras This refers to the Hilbert subspace 
directly described by the BA solution [Therefore, following the studies and notations 



of Ref. [11511 , we can use the a, 7 = pseudoparticles instead of the a, (3 pseudoholes (with 



(3 = ±~) required for the description of the non- LWS's outside the BA solution.] 



The PPT introduced in detail in Ref. [15] involves infinite branches of pseudoparticles 



labelled by the quantum numbers a and 7. Here a = c,s and 7 = 0, 1, 2, 3, ... 00. Fortu- 



nately, for most Hamiltonian eigenstates of physical interest only a small number out of the 
infinite available a, 7 bands have pseudoparticle occupancy. The possible pseudomomentum 
occupancies correspond to the different LWS's of the model. The pseudomomentum values 
q are such that q^ < q < q£~), where the expressions of the pseudo-Brillouin-zone limits, 
q^), are given in Ref. [[TJJ. They are of the form q^ = ±q aa + 0(1/N a ), where for the case 
of a ground state q an is given by 

<? c ,o = 7T ; q ca = n-2k F 7 > , 

q s ,o = k F ^ ; q sn = k F ^-k Fl 7 > . (2) 

The branches c, and s, have been called in previous low-energy studies c and s, 
respectively ||,[18] . They were shown to describe the low-energy excitations of the Hubbard 



chain and to determine the low-energy behavior of its charge and spin transport properties 



18] . (In the limit of low energy the heavy-pseudoparticle branches are empty.) On the 
other hand, description of the LWS's of the model that have a finite-energy gap, u , relatively 
to the ground state involves the heavy pseudoparticle branches c, 7 > and s, 7 > ]H)| |. 

A very useful concept in this theory is that of generalized ground state (GGS). In Ref. 
rj| it was defined as the Hamiltonian eigenstate(s) of lowest energy in the Hilbert subspace 
associated with a given sub-canonical ensemble. The concept of sub-canonical ensemble fol- 
lows from the conservation laws of the a, 7 pseudoparticle numbers, N a ^. Each Hamiltonian 
eigenstate has constant values for these numbers and a sub-canonical ensemble refers to a 
given choice of constant N an numbers. 



On the other hand, in Ref. Jl6[ the concept of GGS was extended to (i) filled 0,7 

-) 



pseudoparticle seas with compact occupations around q = 0, i.e. for qp^y+i < q < QfI j+i 



where the pseudo- Fermi points are given by <?^ 7j+1 = ± ^°' 7 + 0(1/N a ), and (ii) filled a, 7 
pseudoparticle seas with compact occupations for q^ < q < q F a.-y -1 anc ^ f° r Qf. 



(+) 



'a,7,-l — 



q < qij, where the pseudo-Fermi points are given by q^pl^-i — ^[Qa^ — 7r ^°' 7 ] + 0(1/N a ). 
From the studies of Refs. |TB| , [rT| ], it will be shown elsewhere that the creation of one 0,7 



pseudoparticle from the ground state involves, to leading order, a number 27 of electrons. 



Since the currents are two-electron operators, it follows that the creation of single a, 1 
pseudoparticles from the ground state are the most important contributions to the transport 
of charge and spin at finite energies. On the other hand, since the states j^\GS) which define 
our Hilbert space have zero momentum (relatively to the ground state) and the creation from 
the ground state of single a, 1 pseudoparticles of type (ii) is a finite-momentum excitation, 
for simplicity in this paper we restrict our study to GGS's of type (i). Therefore, in order 
to simplify our notation we denote the pseudo- Fermi points qpart+i simply by q F ^ a - These 
are given by qp^~ = ±?f«, 7 + 0(1/N a ) where the pseudo-Fermi momentum [|T5[ 



ixN 

" J, o,7 /q\ 

<?Fa, 7 — —77 , \0) 

1 "a 

appears in several expressions below. Note, however, that the generalization of our results 
to GGS's of type (ii) is straightfoward. We emphasize that Ref. [0 definition of GGS refers 
to the above choice (i) for the c, and s, 7 branches and to the choice (ii) for the c, 7 branch 
with 7 > 0. Therefore, in the case of the c, 7 > pseudoparticles our GGS choice differs 
from the choice of that reference. 

The ground state is a special case of a GGS where there is no 0,7 > heavy- 
pseudoparticle occupancy JT5] and the pseudo-Fermi points (^) are of the form 



q Fcfi = 2k F ; q Fsfi = k Fl ; q Fan = 7 > . (4) 

The PPT consists in expanding the Hamiltonian (||) in the density of excited pseudopar- 
ticles relatively to the initial ground state. This allows us to write Hamiltonian ([!]) in normal 



order relatively to that ground state as [15 



: H := Hq + H Landau , (5) 
where up to second order in the density of excited pseudoparticles, HLandau is of the form 

H Landau = # (1) + # (2) , (6) 

with 



H {1) = E : Na,j{q) ■ , (7) 



= <£?ofa) - e£? («F«,o) , e aa {q) = <M{q) - e£(0) , (8) 

and -f/^ 2 ) is given by 

H {2) = ^E E ;U >7; a' l7 '(g,?') : A\ 7 (?) :: A> Q ,, 7 ,(g') : . (9) 

Here N a)1 (q) = a ~&q,a,7 is the a, 7-pseudoparticle operator number at pseudomomentum g 
and the operators frq jQ7 and & 9 , a , 7 obey the usual anticommuting algebra ]T3| . The Hamilto- 
nian eigenstates decribed by the BA are also eigenstates of N a ^(q) with eigenvalue A ir Q , 7 (g) 
and eigenstates of : N a ^(q) : with eigenvalue 6N a ^(q) = N a>1 (q) — N®^(q), where N® (q) 
is the ground-state pseudomometum distribution These distributions characterize the 



occupancy configurations of the pseudomomenta in the a, 7-pseudoparticle bands. 



The physical meaning of the Hamiltonian terms Hq and HL an d au is explained in Ref. ||15 
These Hamiltonian terms are such that [H , H Landau ] = and H has eigenvalue uj given 
by PJTJ 



uj = 2fi E lN cn + 2/i # EC 1 + 7)^,7 + E < 7 (0)^«,7 > ( 10 ) 

7>0 7>0 a,7>0 

where Af Qi7 are the numbers of a, 7 heavy pseudoparticles created in the transition from the 
ground state to the GGS. 

The set of energies u>o, Eq. (|10|) , play a central role in the theory. This is because for a 
given initial ground state the PPT is associated with a final Hilbert subspace characterized 
by a set of finite N a>1 numbers. The states which span such subspace differ from the initial 
ground state by a small density of pseudoparticles and have small positive (u> — u ) energy. 
The main point of the PPT is that for low-excitation energy, (u — uj ), only the first two 



Hamiltonian terms, Eq. (|J), are relevant |T5|JT6[1 . Therefore, the truncated Hamiltonian 



- (H) describes the physics for energies just above the set of finite-energy values ujq of the 
form flT0|). [This includes u> = 0.] 



Since the conservation of the electron numbers imposes the following sum rules on the 
numbers N an [|15j 

^l = E7^ c , 7 + E( 1 + 7)^, 7 , (11) 

7>0 7 

and 

N = N cfi + 2j2lNc^, (12) 

7>0 

the creation of heavy pseudoparticles from the ground state at constant electron numbers 



requires the annihilation of a, pseudoparticles. It follows from Eqs. ( p)) and (|12|) that the 



changes AN at0 associated with a corresponding ground-state - GGS transition read 

AiV s , = - E 7^ c , 7 - EC 1 + 7)iV s , 7 , (13) 

7>0 7>0 

and 

AiV Ci o = -2^ 7 iVc,7. (14) 

7>0 

For instance, the creation of one c, 7 heavy pseudoparticle from the ground state requires the 
annihilation of a number 27 of c, pseudoparticles and of a number 7 of s, pseudoparticles, 
whereas the creation of one s, 7 pseudoparticle involves the annihilation of a number 1 + 7 
of s, pseudoparticles and conserves the number of c, pseudoparticles. 

Although, following Eq. (|10D, ujq can be large, we emphasize that the PPT is always 
a low (ui — uj ) energy theory. This is because within the PPT the densities of removed 
a,0 pseudoparticles, —AN at0 /N a , of added 0,7 > heavy pseudoparticles, AN a>J /N a = 
N a ^/N a , and of their pseudoparticle - pseudohole processes are always kept small. Moreover, 
for each set of finite iV Qi7> o numbers there is one PPT and one value of energy (|T0|). 

The energy bands e° a Jq) and the /-functions /o i7;a ',y(g,g') are given, respectively, by 
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4,o(l) = 2t™*K {0) iQ) + 2t J Q *$ c ,o ;Cl o (k, K (0) (g)) sinfc , (15) 



S 



< 7 (g) = -1 U + 4tRe\fl-u 2 [I&%)-i7] 2 + 2t f Q rffc$ c , 0;c , 7 R®(q)) sin fc , (16) 
e° 7 (g) = 2t / Q ^$ Ci0 ;s, 7 (A;, i$>(g)) sin fc , (17) 

and 

^ fa, r ,a',y(q, q') = w a , 7 (g)$ a , 7;a ', 7 '(g, g') + w a / i7 /(g')$ a / i7 / ;a)7 (g / , q) 

oo 

+ ^2 ^ ^ 0(N a " ^")v a " ^"^ a " ^" :a ^(jqFa" ^" , q)®a",j";a',y(jqFa",y, gQ • (18) 
j'=±l a" 7" 

Here = 1 for x > and #(2) = otherwise is the usual theta function and v an = 
u a , 7 (gFa,7) * s ^he velocity at the pseudo-Fermi point. The phase-shift functions $a i7;a ' )7 ' 
and the phase shifts $a, 7 ;a', 7 ' are defined in Ref. [0. ^^^^/(g, g') measures the shift 



in the phase of the a', 7' pseudoparticle of pseudomomentum q' due to a zero-momentum 
forward-scattering collision with the a, 7 pseudoparticle of pseudomomentum q. It is useful 
to introduce the function W°(q) such that W = K, R cr/ , R sr/ . (Here, 7 = 1, 2, . . . , 00 for 
R crt and 7 = 0, 1, 2, . . . , 00 for R s>7 .) It represents any of the three ground-state rapidity 
functions K^°'(q), R^(q), and R^(q), whereas the functional W(q) represents any of the 
three general functional rapidity functions K(q), i? Cj7 (g), and i? S)7 (g) |0|. These functional 
are obtained from Eqs. (|A7|) , (|A8|) , and (|A9|) with = 0. (In that Appendix we solve the 
BA equations with twist angles.) The ground-state functions W°(q) are obtained by taking 
the particular choice N a ^(q) = iV° (g) in the latter equations. It is useful to introduce the 
pseudo-Fermi rapidity parameters [15[ 



Q = K(°\q Fc>0 ); r c , = ^ ; r a , 7 = R®(q Fa , 7 ) , (19) 

where Q appears in the integrals of Eqs. ( |T5| ) - (|T7|) and in the last expression r an refers to 
all a, 7 branches except c, 0. 



III. CHARGE AND SPIN CURRENTS: SOLUTION OF THE BA EQUATIONS 



Within linear response theory the charge and spin currents of the ID Hubbard model 
can be computed by performing a spin-dependent Peierls-phase substitution in the hopping 
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integral of Hamiltonian (0), t -> te M ^l Na |||27 . 

It has been possible to solve the Hamiltonian ([!]) with the additional hopping phase 



3±i4>a/N a 



by means of the coordinate BA both with twisted and toroidal boundary conditions, 



both approaches giving essentially the same results [pTjpg]] . One obtains the energy spectrum 
of the model parameterized by a set a numbers {kj,As} which are solution of the BA 
interaction equations given by 



N 



i h tt sin(fcj) - A<5 + iU/A 
£ sin(%) - A* - iU/4 



N), 



(20) 



and 



" sin(fcj) -A s + iU/A 
W sin(A; J ) - A 5 - iC//4 ~ & 



n 



A^g — Aj - i?7/2 



(5 



(21) 



However, previous studies of the CT 7^ problem p7|j28|1 have only considered the real 
BA rapidities of Eqs. ( p0|) and ( pT|) which refer to low energy. Here we follow the same steps 
as Takahashi ||J] for the CT = interaction Eqs. fl2"0| ) and fl2~T|) and consider both real and 
complex rapidities. We then arrive to the following CT 7^ equations which refer to the real 
part of these rapidities 



k,N a = 2711" 



N c . 



EE 2to 

7 >0 j'=l 



EE 2 tan 

7 j'=i 

_ 1 / sin(/c J )/w — _R, 



_ x / sm(kj)/u — R 



(7 + 1) 



c,7J 



7 



(22) 



_ x / sm(kji)/u — -R c ,7,j 



2N a Re sin _1 ([ii; C)7)j + ij\u) = 2tt/ j c ' 7 + 7(0 T + ; ) - ^ 2 tan" 



^.y 



+ E E ®7:7'(-^c,7,j RctYd' 

7 '>o i'=i 



and 



sin(A;')/-u N 



E^tan-M ^ :T J// " ) +(l + 7)(0T-0l) 



(1 + 7) 



N. 



(23) 



2 ^* 7 + E E ©7+i,7' + i(^,7 J - Rsn',f) ■ (24) 



y j'=i 
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The functions 6 7i y(x) [and 7+ i j7 /+i(#)] of Eqs. (p2|), fl23|), and fl24|) are defined in Ref. |L5 . 



The following definitions for the real part of the rapidities, A™ +l /u = R s ,j,j (with n + 1 = 7 
and a = j), A'j l /u = R c ,-y,j (with n = 7 and a = j), and 7 = 1, 2, . . . , 00 for the iV ( 



sums 

and 7 = 0, 1, 2, . . . , 00 for the iV Si7 sums, allows us to recover Takahashi's formulae for = 
U|. Here and often below we use the notation c = c, 0, which allows the c, 7 sums to run 
over 1,2,3,..., 00. Whether we are using this notation or the previous one will be obvious 
from the context. 

The important numbers IJ, Ij' 7 , and J|' 7 which appear in going from Eqs. (|20|) and 
(pl|) to Eqs. (0), (p3|) , and ( p4|) are the quantum numbers which describe the Hamiltonian 
eigenstates. Depending on the parity of the numbers [H 7 =o N s ^ + J2j=i N c>7 \, [N a — N+N cn ], 
and [N — N Sjl ], respectively, they are consecutive integers or half-odd integers |15| . All the 
LWS's of the model are described by the different occupancies of these quantum numbers. 
For example, the ground state is described by a compact symmetric occupancy around the 
origin of the numbers 7? and Ij' , and by zero occupancy for the numbers Ij' 1 and J|' 7>0 
|T5| . It is convenient to describe the eigenstates of the model in terms of pseudomomentum 
{qj n = 27r/j*' 7 /N a } distributions, where If = Ij. 

The energy and momentum eigenvalues are given by |T5| 



N c N c . 



E = -2tJ2 cos(%) + 4tRe^Jl - u 2 [R c ^j - i 7 ] 2 

j=l 7>0 j=l 

+ N a (U/A -fi) + N(jji - U/2) - ^H{N^ - N{), (25) 



and 



3=1 7 j=l 7 >0i=l 



vr ^ K n , (26) 

7>0 



respectively. We emphasize that the rapidity dependence on </> CT is defined by Eqs. 
and determines the energy-functional ( p5j) dependence on <p a . The corresponding <p a = 
expressions recover the rapidity and energy expressions of Ref |15| . 



In the limit of a large system (N a — > 00, N/N a fixed) we can develop a generalization 
of the low-energy pseudoparticle-Landau-liquid description of the Hubbard model [pOH 22 
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and of its operational representation |18|]. This generalization refers to energies just above 



the set of energies co> , Eq. (|10l), where the Hamiltonian (|5[) describes the quantum-problem 
physics. Note that the choice uq = 0, which refers to N a>1 = for 7 > 0, recovers the usual 



low-energy theory of Refs. p8| . [20| -pl| . On the other hand, when luq > 0, in addition to finite 
occupancy of the usual c, = c and s, = s pseudoparticle bands, there is finite occupancy 
for some of the branches of the heavy c, 7 and s, 7 pseudoparticles [|15[ . 

In the above thermodynamic limit the rapidity real parts kj = kj(qj), -R s , 7 j = R S}J} j(qj), 
and R c ,-y,j = Re,-y,j(qj) proliferate on the real axis. As in Refs. pOH22|l, equations (p2|) , 
(p3j), and (pif) can be rewritten as integral equations with an explicit dependence on the 
pseudomomentum distribution functions N a ^(q). These are Eqs. (|A7|), ( |A8|) , and ( |A9|) 



of Appendix [A] which refer to the case CT 7^ 0. In that Appendix we derive ground-state 
normal-ordered expressions for the rapidities and charge and spin currents. 

The combination of Eqs. (p5|) , ( |A7|) , ( |A8|) , and ( |A~^ ) allows the evaluation of several 



interesting transport quantities. This includes the charge and spin currents and the charge 
and spin pseudoparticle transport masses. The charge and spin current operators (with 
C = p for charge, and ( = a z for spin) are for the ID Hubbard model given by 

N a 

f = -eit J2 H(c]o9+io - c] +la c ja ) , (27) 

a j=l 

and 

N a 

Y' = -(l/2)«E&(^+to " ■ (28) 

cr j = l 

Importantly, the discrete nature of the model implies that the commutators of the Hamil- 
tonian ([!]) and of the current operators j^, Eqs. fl2~7|) and (f2~8|), are non zero. It follows 
that the BA wave function does not diagonalizes simultaneously the Hamiltonian (JI|) and 
the current operators ( |27|) and (p8|). Since the BA solution alone only provides the diagonal 
part in the Hamiltonian-eigenstate basis of the physical operators , we can only evaluate 
expressions for the diagonal part of the currents which provide the mean values of the charge 
and spin currents. These refer to all LWS's and are important quantities for they allow us 
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to compute the transport masses of the charge and spin carriers of the system. In addition, 
our formalism defines the charge and spin carriers. These are found to be the c (i.e., c, 0) 
and c, 7 pseudoparticles for charge, and the c and s, 7 pseudoparticles for spin. This follows 
from Eqs. ( |A7| )-(|M3|) and also from the Boltzmann transport analysis of Sec. [V]. 

We emphasize that combining the generalized pseudoparticle representation |TJ| with 



a low-energy (u — ujq) conformal- field theory []lq1 , leads to finite-energy current - current 



correlation function expressions which are determined by the non-diagonal terms (in the 
Hamiltonian-eigenstate basis) of the current operators. This is a generalization of the low- 
energy correlation-function studies of Refs. |j,|23|. However, these expressions cannot be 



derived within the BA solution alone. Therefore, these studies go beyond the scope of the 
present paper and here we consider the diagonal part of the charge and spin current operators 
only. 

The mean value of the current operator in a given LWS, |m), is given by 

d{E m /N a 



(m\ j^\m) = — 



d(<P/N a 

where E m is the energy of the Hamiltonian eigenstate \m) and ]27] 



(29) 

0=0 



= T = , ( = p, 

= T = -0 i , ( = cr z . (30) 

In our basis the LWS's are simply obtained by considering all the possible occupation 
distributions of the pseudomomenta qj = 2wl c * a /N a . Therefore, it is convenient to describe 
the matrix elements (m\j^\m) in terms of the pseudomomentum occupation distributions 
N a ^(q). This leads to a functional form for the current mean values. The computation of 
(m\j^\m) involves the expansion of Eqs. ( p5]) and (|A7|) - (|A9| ) up to first order in the flux 0. 



Writing Eq. (p5|) in the limit of iV a — > 00, expanding it up to first order in the flux 0, and 
using Eq. (|29|) we obtain 



(m\j c \m) = -2t— r dqN c {q)K' t> {q)sm{K{q)) 

2lT J-g r _ 



2lT J-q c 
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E 4 V" n dqN ca (q)Re 

Z7T J —q c ,~t 



(?), 



(31) 



7>o "~ yc ^ \Jl- u 2 [R cn (q) - ryp 

where the important functions W^(q) (with W = K, R S}7 , and R c ,y) are the derivatives of 
the rapidity functions defined by Eqs. ( |X7| ) - in order to the flux <fi at = 0. They 
obey a set of integral equations obtained from differentiation of Eqs. ( 1^7 ) - ( |A9| ). 

It is convenient to write (m| j^m) in normal order relatively to the ground state. To 
achieve this goal we expand all the rapidities W(q) and the functions W^(q) as 



W(q) = W\q) + W\q) + ... , 
W+(q) = W°>*(q) + W W (q) + 



(32) 
(33) 



respectively. In these equations the functions W°(q) and W 0, ^(q) are both referred to the 
ground state, and the functions W 1 (q) and W 1, ^(q) are first-order functionals of the devia- 
tions 8N a ^(q). In Appendix |A] we show that the above expansions lead to a ground-state 
normal-ordered representation. To first order in the deviations the normal-ordered expres- 



sion for the matrix element (|3~TD simply reads 

(m\j C \m) = ^2^2 dq5N a „(q)ji„(q) 



(34) 



a 7 



where the elementary-current spectrum j^Aq) is given by 



Ja,7 



(35) 



a' 7' 



Here 



-F«,7; a ',7'('?) - — X! i/<*,7;a' ,7' (<?,:/' <?Fo 



(36) 



and C£ 7 are the coupling constants of the pseudoparticles to charge and spin given by 



c l 7 = 8 a ,c$j,o + 



a,7 i 



(37) 



where 



K£„ = <L C 2 T ; 



a,7 



/C' 



a, 7 



-5 QiS 2(l + 7 ) 



(38) 
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As in a Fermi liquid |3l|,32|, the expressions of the elementary currents (|35|) involve the 
velocities v an (q) and the interactions [or /-functions] fa^a^yiq, ?')• However, the pseu- 
doparticle coupling contants to charge and spin, Eqs. ( pTf ) and fl38|), are very different from 
the corresponding couplings of the Fermi-liquid quasiparticles. We emphasize that at low 
energy Eq. (^) recovers the expression already obtained in Ref. || which only contains 
the c = c, and s = s, elementary currents. The coupling constants (|37|)- (|38|) play an 
important role in the description of charge and spin transport and are a generalization for 
7 > of the couplings introduced in Ref. ||. They define the a, 7 pseudoparticles as charge 
and spin carriers. We emphasize that when C^ 7 = the corresponding a, 7 pseudoparticles 
do not couple to ( (i.e. charge or spin). Therefore, for 7 > the 0,7 and s,7 pseudopar- 
ticles do not couple to spin and charge, respectively. This is related to the charge and spin 
separation of one-dimensional quantum liquids which in the case of the present model was 
studied in Refs. fL8p4l . Importantly, when C£ = the 0,7 pseudoparticle - pseudohole 



processes do not contribute to the ( correlation functions. It will be shown elsewhere that 
this provides a powerful selection rule which implies that some of the terms obtained from 
the small (u — Uq) conformal-field theory |]TB|| for the correlation functions vanish. 



In contrast to the general current expression ([H]) , expression (|3^) is only valid for Hamil- 
tonian eigenstates which differ from the ground-state pseudoparticle occupancy by a small 
density of psudoparticles. This is because in expression fl34]) we are only considering the 
first-order deviation term. 

The velocity term of current-spectrum expression (|35|) is what we would expect for a non 
interacting gas of pseudoparticles and the second term takes account for the dragging effect 
on a single pseudoparticle due to its interactions with the other pseudoparticles. (This 
is similar to the Fermi-liquid quasiparticle elementary currents We remind that 

Eq. ( p5|) is valid for finite energies ui just above the energy ui Q corresponding to the suitable set 
of finite N an numbers. These numbers characterize the state \m). Therefore, the sum over 7 
is in Eq. Q55D restricted to the a, 7 bands that have non-zero occupancy of pseudoparticles, as 
is imposed by the presence of the step-function. Within the PPT, the deviation second-order 
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pseudoparticle energy expansion corresponds to the deviation first-order current expansion 
([53]) which refers to small positive values (u — u>o) of the excitation energy. In contrast to 
Fermi liquid theory, our PPT is valid for finite energies [just above the energy values c^o, Eq. 
(^0|)1 because (i) there is only forward scattering among the pseudoparticles at all energy 
scales and (ii) at small {oj — uj q ) energy values only two-pseudoparticle forward-scattering 
interactions are relevant. (In a Fermi liquid this is only true for ujq = and uo — > [pl] , |32|| .) 
We emphasize that the current expression (|3lD includes all orders of scattering and, therefore, 
applies to all energies without restrictions. Elsewhere it will be shown to be useful in the 
study of charge and spin transport at finite temperatures. 



IV. PSEUDOPARTICLE TRANSPORT AND STATIC MASSES 

In reference || the charge and spin transport masses of the c, and s, pseudoparticles 
were defined and were shown to play an important role in the transport of charge and 
spin. For instance, they were shown to fully determine the charge and spin stiffnesses 
P JT0|JT3||27|j29|j30[| . Here we generalize the mass definitions of Ref. to 7 > and define 



the charge and spin transport masses, , as 



where j£ = i£ 7 (gi? a , 7 ). They contain important physical information. As in a Fermi liquid 
31,32fl, the ratio 



r l 7 = m i 7 / m «,7' ( 40 ) 

of the transport mass over the static mass provides a measure of the correlations importance 
in transport. Similarly to the 7 = case 0, the latter is in general defined as 

< 7 = ^ • (41) 



In Appendix [B] we define the mass ([II) in terms of suitable functions and find some limiting 
expressions. 
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It can be shown from the transport- and static-mass expressions that the ratio m£ 7 /m* 
involves the Landau parameters 

F a,r,a',l' = F l, r ,a.'rt'( < lF<*,'l) ' i = °> 1 ' ( 42 ) 

with , ,(q) given by Eq. (|36|). These parameters can be written as follows 

F a, T ,a',Y = -^a,a'^,j'V a>1 + Y @{ N a>> ,Y>)v a " ,-y" \C a " ^" - an C a " rf " - a ' n '} , (43) 

a". 7" 

where the quantities C a , r ,a',y are given by 

C 7 ;a',7' = <W'&y,7' + ? ® a, r ,a.' (<lFa,j, ^Fa'.f) ■ (44) 

We find for the ratios Tn^/m^ the following expressions 

7 = 0, (45) 



!,7' a <7 
m «,0 u a,0 



Q,0 



Ca, 7 (So!',a" ^Q',O t '«",o(a»0;a,0^",0;a',i 



and 



777C 1 

-^ = - — c < — 1 ' ^ >0 - ( 46 ) 

m a,7 Ca, 7 (C a ,7 + Y,a' ^a',O^Q,7;a',o) 



In the Table analytical limiting values for the mass ratios of form (]4"0 ) are listed. Obviously, 
since for 7 > the c, 7 and s, 7 pseudoparticles do not couple to spin and charge, respectively, 
the ratios m a c z ^j m* and m p s / m* are infinite. 



As was referred previously, it can be shown from the results of Refs. U1~5}-|1~71 that the cre- 
ation of one a, 7 pseudoparticle from the ground state is a finite-energy excitation which, to 
leading order, involves a number 27 of electrons. Therefore, and since the current operators 
are of two-electron character and couple to charge and spin according to the values of the 



constants (|37D and (p8|), at finite energies the c, 1 and s, 1 heavy pseudoparticles play the 
major role in charge and spin transport, respectively. On the other hand, the a, 7 > 1 heavy 
pseudoparticles contribute very little to charge and spin transport. It follows that in the 
present section we restrict our study to the ratios ([40]) involving 7 = 1 heavy pseudoparticles. 
We consider the ratios m p cl /m* l and rrf s y/m* v (Note that rn°\/rn* cl = m p sl /m* sl = 00.) 
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We also consider the case of the c, charge-mass ratio which is closely related to the charge 
stiffness studied in detail in Ref. 0. In Figs. [THT^ these ratios are plotted as functions of the 
onsite repulsion U in units of t, electronic density n, and magnetic field h = H/H c . Note 
that the ratios of the figures are smaller than one. Moreover, the a, 1 mass ratios never 
achieve the value 1 whereas the a, mass ratios tend to one in some limits because of the 
generalized adiabatic principle of Ref. §| . 

Combined analysis of Figs. [I] and @ reveals that the charge-mass ratio for the c, 1 pseu- 
doparticle is, for large U, fairly independent both of the band filling n and magnetic field h. 
It is a decreasing function of U and as a function of the density, n, goes through a maximum 
for a density which is a decreasing function of U. Moreover, figures |3] and f| show that this 
ratio is a decreasing function of the magnetic field. 

In contrast, Figs. |5| - [8| reveal that the charge- mass ratio for the c, pseudoparticle is an 
increasing function of U and of h and as a function of the density, n, goes through a minimum 
for a density which is a decreasing function of U. Note that from Fig. |5| the evolution of 
the c, pseudoparticles to free spinless fermions as U increases is clear. This is signaled 
by the ratio going to one as U — > oo. This behavior follows from the generalized adiabatic 
principle of Ref. [[| and agrees with the well known decoupling of the BA wave function in 
free spinless fermions (in the low-energy sector [flOU ) and localized antiferromagnetic spins 
||41|| . Figures |7| and || also reveal that in the fully-polarized ferromagnetic limit, h — *■ 1, 
the ratio goes to one. This mass-ratio behavior also follows from the generalized adiabatic 
principle || and confirms that in that limit the onsite Coulomb interations play no role in 
charge transport (they are froozen by the Pauli principle) . 

Note that in the large-?/ Figs. - (c) and || - (c) the ratios m p cl /m* cl and Tn^ /m* , 
respectively, are almost symmetric around the density n = 0.5. This implies that for large 
U the charge transport properties show similarities in the cases of vanishing densities and 
of densities closed to one. 

Figure |S| shows that the spin-mass ratio of the s, 1 pseudoparticles is a decreasing function 
of U but that it depends little on U for U > 6. For large U this ratio almost does not depend 
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on the density n, as revealed by Fig. |T0| — (c). Figures |TTj show that, in general, it is an 
increasing function of n but that for h — *■ 1 it has a maximum for an intermediate density. 



In figures |11| and [12] this spin-mass ratio is plotted as a function of h. It is a decreasing 
function of h 

The transport masses are very sensitive to the effects of electronic correlations, as for 
instance to the metal-insulator transition which occurs at zero temperature when n — > 1 |§ . 
As a direct result of this transition, m p C Q — > oo as n — > 1, as was shown and discussed in 
Ref. |§. Moreover, the zero-temperature charge and spin stiffnesses, D^, P. [lO| . [T3| , ^7| , |29| , |30 
defined as 



dC _ l d 2 (E /N a ) 

2 d^/Naf 



(47) 

</>=0 



where <fi is defined for charge ( = p and spin ( = <j z in Eq. (|30|), are such that 2ttD^ = 
J2 a lFa,o/ m a,o • F° r charge, m^ = oo, and the latter expression reads 2irD p = qFc,o/ m c,Q 
and is such that D p — > as n — > 1, satisfying Kohn criterion These quantities can 



be computed within the pseudoparticle formalism by direct evaluation of Eq. (£7). They 



can also be obtained by combining a pseudoparticle Boltzmann transport description with 
linear response theory, as in Ref. [[| . In order to confirm the validity and correctness of our 
formalism, in Appendix |C| we have recovered the charge and spin stiffness expressions (135) 
- (137) of Ref. by direct use of Eq. (0). 

Equations (|37|) and (p8|) show that the Hubbard-chain charge carriers are the c, 7 pseu- 
doparticles. In contrast to the zero-temperature limit where the c, pseudoparticles fully 
determine the charge stiffness, we expect that the c, 7 heavy pseudoparticles play an impor- 
tant role in the charge-transport properties at finite temperatures 0,0. Moreover, else- 
where it will be shown that the limiting behavior of the s, 7 and c, 7 heavy-pseudoparticle 
bands as H — > and n — > 1, respectively, will have effects on the charge- and spin-transport 
properties at finite temperatures. In order to obtain some information on that behavior, 
it is useful to consider limiting values for the quantities whose general expressions we have 
introduced in previous sections. 
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In Appendix [D] we present simpler equations to define the pseudoparticle bands and 
phase shifts in the limit of zero magnetic field. These results show that for H — > and 
7 > the bands e° 7 (g) collapse to a point. This is because both the bandwidth [see Eq. 
([EXE])] and the momentum pseudo-Brillouin zone width [see Eq. (Q)] go to zero as H — > 0. 
This behavior is also present in the Heisenberg chain and, therefore, in that model the 



triplet and singlet exitations are degenerated at zero magnetic field ||39|| . This also holds 
true for the Hubbard chain at H = and in the limit U ^> t, where the BA wave function 
factorizes in a spinless-fermion Slater determinant and in the BA wave function for the ID 
anti-ferromagnetic Heisenberg chain. On the other hand, in the limit n — ► 1 the bands 
e° 7 (g) (for 7 > 0) collapse to a point also because both the bandwidth and the momentum 
pseudo-Brillouin zone width [see Eq. (0)] go to zero in that limit. 

V. KINETIC EQUATIONS FOR THE PSEUDOPARTICLES 

In the previous sections the quantum-liquid physics for energies just above the ujq val- 
ues, Eq. fllUD , was described in terms of homogeneous pseudoparticle distributions. The 
pseudoparticles experience only zero-momentum forward-scattering interactions at all en- 
ergy scales. This property is absent in Fermi-liquid theory where it holds true only at low 
excitation energy when the quasiparticles are well defined quantum objects |]3l|,|32|,|37| . This 



unconventional character of integrable models [38] allows us to extend the use of the kinetic 



equations to energy scales just above the ojq energy values, Eq. fllOD , and not only to low 
energies @. The results presented in this section are a generalization of the kinetic-equation 
low-energy study presented in Ref. 0. 

In the final Hilbert subspace of energy u relative to the initial ground state the Hubbard 



model can be mapped onto a continuum field theory of small energy (lo — ujq) [|16|]. The 
time coordinate t of such theory is the Fourier transform of the small energy (u — loq) which 
corresponds to a finite energy uj in the original Hubbard model. The validity of this approach 
is confirmed by the fact that it fully reproduces the rigorous results of Section |TT[ 
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Let us consider excitations described by space and time dependent pseudoparticle dis- 
tribution functions, N a ^(q,x,t), given by 

N aa (q, x, t) = iV° 7 (g) + 5N Q ^(q, x, t) , (48) 

where N® (q) is the ground-state distribution. It follows from the PPT introduced in Ref. 



1 15] and discussed in the previous sections that the single-pseudoparticle local energy is 



given, to first order in the deviations SN ajl (q,x,t), by 

£ a , 7 (g, x, t) = e a , 7 (g) + — dq5N a , a ,(q , x, t)f a ^ a , : y(q, q) . (49) 



Let ^ represent the total charge, ( = p, or spin, £ = cr z . It follows from the relations 
(P"T|) and (IT) involving the pseudoparticle and electron numbers that AS depends linearly 
on the pseudoparticle deviation numbers. Thus, in the case of inhomogeneous excitations 
described by Eq. fl48|) the corresponding expectation value at point x and time t, (A''(x,t)), 
can be written as 

(A<(x,t)) = (A% + ^ £ r" 1 ' dq'5N a ,, Y (q',x,t)Ci, tY x of, (50) 

where a p = —e and a" z = 1/2. 

In this "semi-classical" approach the response to a scalar field, V*(x,t), is proportional 
to the conserved quantity A?. As for low energy ||, in the presence of the inhomogeneous 
potential the force ^(x, t) a ^ that acts upon the a, 7 pseudoparticle is given by J 7 ^ (x, t) = 
— [dV>(x, t)/dx]C^ x aS. It follows that the deviations 6N a ^(q,x,t) are determined by the 
solution of a system of kinetic equations (one equation for each occupied a, 7 branch) which 
reads 

Q _ dNg tl (q,x,t) dN aa {q,x,t) de aa {q,x,t) _ dN a „{q,x,t) de a „(q,x,t) 
dt dx dq dq dx 

dN aa {q,x,t) dVt{x,t) nC c 

dq dx~ C ^ X a ■ (51) 

Introducing Eq. (fPs|) in Eq. (|5Tj) , expanding to first order in the deviations 5N a ^(q,x,t), 
and using Eq. fl49D we obtain the following set of linearized kinetic equations 
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d5N a „(q,x,t) . ^ d6N a „(q,x,t) 

dt + v <vrW dx 

d5N an (q,x,t) f dV<(x,t) 



dq \ dx 



a, 7 



1 /-W A , d6N a r )Y (q',x,t) ,1 

The conservation law for (^(x, i)) leads in one dimension to 

d{A<(x,t)) (J<(x,t)) 

—di + ^^ = °' (53) 

where (A^(x,t)) is given by Eq. (|57JD and (J' ( '(x,t)) is the associate current. Multiplying 
Eq. ( |52|) by C£ x a^, summing over a and 7, and integrating over g we find for V^(x, t) = 
and by comparing the result with Eq. fl53|) that the current spectrum (g) is given by 
cp times expression (|35D . (This expression has been derived from the solution of the BA 
equations with = 1.) 

This agreement confirms the validity of the above low- (a; — uj ) continuum- field theory. 
The unusual spectral properties associated with the zero-momentum forward-scattering char- 
acter of the pseudoparticle interactions follow from the integrability of the Hubbard chain 
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VI. CONCLUDING REMARKS 



In this paper we have generalized the finite-energy PPT [T5|] to the case of the Hubbard 



chain with a spin dependent Peierls substitution. This has allowed the evaluation of the 
charge and spin currents in terms of the elementary currents of the charge and spin carriers. 
We have shown that at all energy scales these carriers are the a, 7 pseudoparticles of the 
PPT. We have evaluated their couplings to charge and spin and introduced the associate 
charge and spin transport masses. Our results are also a generalization for finite energies 
of the low-energy studies of Ref . , our charge and spin current expressions recovering the 
expressions already obtained in that reference in the limit of low energy. 
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The obtained heavy-pseudoparticle transport masses are important quantities. They 
are believed to control the charge and spin stiffnesses at finite temperatures. Moreover, the 
pseudoparticle couplings to charge and spin obtained in the present paper provide important 
selection rules concerning the ground-state transitions which contribute to the finite-energy 
charge - charge and spin - spin correlation functions. Expressions for these functions can be 
derived by combining the BA solution with a low-energy (u — ujq) generalized conformal-field 
theory ||16|| . The above selection rules will be shown elsewhere to provide important informa- 



tion on the finite-energy correlation functions which cannot be extracted from conformal-field 
theory alone. 

Finally, the CT -dependent charge and spin current expressions of general form ([31]) will 
be used elsewhere to find out whether the half-filling Hubbard model is or is not an insulator 



at all temperatures 34,35 
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APPENDIX A: NORMAL-ORDERED SOLUTION OF THE BA EQUATIONS 

WITH THE FLUX <j> 

In this appendix we derive the normal-ordered BA equations required for the evaluation 
of Eqs. (0) and (j35|). Writing W*{q) from Eq. © as 

W%) = , (Al) 

where L equals L Cj0 , L Cj7 , or L Sj7 when W equals K, R Sjl , or R Ctl , respectively, we find that 
W 1,( ^(q) obeys the following equality 
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W»(q) = ^L^(q) + ^°'*(?) . (A2) 



Introducing the above equation in Eq. (plD and writing the distributions functions N a ^(q) 
as N® (q) + SN a%1 (q), we can expand J = (m\j^\m) in terms of the pseudomomentum 
deviations as 



J = J + J + J ... , (A3) 

where the first term, J 1 , of the current normal-ordered expansion ( |A3|) can after some algebra 
be written as 

rl rl 0+ V- L °'oO'^) L c,o(igFc) . . . 

J =J °" 2t J -S, SmW) + 

7>o j=±i Jl — u 2 [ir ca — ij] 2 27rp Cj7 (r Cj7 ) 

_ ^ p dq ^M siniK^iq))^) , (A4) 
J-g Fc dq 

where the functions 2iTp C) o(k) and 2 r Kp olil {r) were defined in Ref. []15 and 
J-qc dq 

+ ^ReAt dq5N^(q)—= ' r — — -f L ^(9) . (A5) 

7>o V 1 - u2 [ i ^'r(9)- i T] 2 g 

The function £ 1,<?i (g) is defined as 



C^(q)=L^(q)-W\q)^l 



(A6) 



In order to obtain the integral equations for L°'^(q) and C 1 ^^) (with C = £ Cj0 , Ay-p and 
C s>1 ), we start from the continuum limit of Eqs. (p2]), (p3|), and fl24|) which reads 



_ E 1 /««' 2tan -. f ^i!^M | , (A7) 
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2Re sin- 1 [( J R Cj7 ( g ) - ij)u) =q + 7 (0 T + (f>i)/N a - 

1 f qc j„iATf„'\n f sin(if(g '))/u - R C)J (q) ' 



- 7; dq'N c (q')2tan 

Z7T J- 0r 



7'>0 " g c, 7 ' 

and 



-<fc V 7 / 

+ E ^ ^'^ c ,y(^)07,7'(^c, 7 (g) - R c ,Aq') » ( A8 ) 



, = ( 7 + l) lh - M ,N. + ± jT ') 2 tan- =^»Zg ) 

- ttq'N„,(q')e^ w (R,, y (q) - «„-(?')) • (A9) 

, ATI J-q s ~ 



It is convenient to write the function ,(x), defined by Eq. (B7) of Ref. fT5fl ) as follows 

2bT 

= ^T+[x/l} 2 



e!vW = E7^' (Aio) 



We emphasize that comparision term by term of expression (B7) of Ref. |T5[ with expression 



( ]A10|) fully defines the coefficients b]' 1 and the corresponding set of integer numbers I. 

Following equation fl30|), we have 0j = 4>± for a charge-probe current and <f>^ = — <fi^ for a 
spin probe. With the above equations written in terms of 0f and <fii, Eq. (|A4| ) provides both 
the charge and spin currents. In what follows, we introduce in the functions L^(q) the index 
( = p,a z to label the equations for either the charge or the spin current, respectively. We 
start by expanding Eqs. fl^?D, fl^5|), and ( |A9| ) up to first order in 0. This procedure reveals 
that the functions L^^(q) obey the following integral equations 

Lc > o[q) - Cc '° , (Y + 1)tt i- 9s , y 5 1 + r HK( g ))/u-R^ ig ') v , dq > L s,A a ) 

+ v — f c ' y d g / Nc ' Yiq,) dR cAq') L Mu (al) (An) 

+ h *i q 1 + pmtz^Sl ? dq> {q } ' (AU) 

+ VV - T C ' V do' N ^'^ d^Ml L U (q/) (A12) 

hi*"-** l+[ ^i3t^l ? dq' L ^ {qh (M2) 

and 
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1 f 9 ° Nc(q') dK{q>) _ 1t ^^ tU( ,. 



1 P-y N sn ,{q')b] +l >y +l dR sa ,{q>) TU( 



- ^ T ^^^ -^r 1 ^ ' (A13) 

where the coupling constants C£ are defined by Eqs. (p7|) and ([381). We again write the 
distributions functions N ar/ (q) of Eqs. flATTD , ([AT2D , and ( |AT3|) as iV° )7 (?) + 5N an (q). This 
allows us to obtain integral equations for L '^(g) and C 1,( ^(q) (we remark that the functions 
C l,, P(q) are the same both for £ = p, a z ). It is then straighforward to find the integral 
equations obeyed by L 0,l ^^(q) and show that L 0,( ^^(q) can be simply expressed in terms of 
linear combinations of phase shifts. The final result is 

L^{q) = Ci n + E E &(N«rf)(*rf*^M(q,3q F ^) . (A14) 

a',7' i=±l 

The integral equations obeyed by C 1,(l> (q) are related to the integral equations obeyed by 
£ 1 ' < ^(r), where r equals sm(K^(q))/u, R^(q), and R^(q) for C = C C:0 , £ Cj7 , and £ S)7 , 
respectively. The functions C x, ^{r) obey the following integral equations 



4tf (r) = ^'°(r) + i P° ^ 5° (r \ 2 > ( A15 ) 



and 



' 7r( 7 + l) M 7- rc l + (^) 2 



L -j / rfr i , rr-r'12 ' ( AlT ) 



where the free terms £^'o' ( r )) Avy ,0 ( r ); an d £5,7'° ( r ) are > respectively, given by 
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£-7-7 



-f 7T(Y + 1) J- q . 



la.y' 



dq'5N S! y{q')- 



r r - R i°!/(9') n 



dq' 



V 7 27rp C) y(r Cj y) (i + [!^]2 ) 



s,7' 



E 



(Y + l)vr^i 27rp s ,y(r S)7 (i + f^f]*) 



(A18) 



7 



and 



7U7T 2vrp c ,o(Q) (1 + p*^]*) C ° S W) 

- E W*y) E ^ E x 2npc ArcY) (1 + R ^ ]2) . 



£iAO (r) = 1 /* ^ cos( K^q')) dKi0) ^ 

^ [ ] nu(rf + l)J- qc q c[q) i | [ g ^MZHzr ]2 C ° SlA WJJ dg/ 

~VttZ /-,..„, 8,7 W \ r r-a?lW„ dq' 



7' J 



1 v jLlotiQFc) L°J(jq Fc ) 

27rp Ci0 (Q) (1 + [1^°]2) C ° SWj 

- E E ^ E x 27rp Si y(r Si y) (l + [^] 2 ) • (A20) 



Introducing the functions C^i^q) obtained from Eqs. (|A15|) , (|A16|) , and ( |A17|) , in Eq. 
and keeping terms only up to second order in the density of heavy pseudoparticles, we obtain 
Eq. (gT]) with ii n (q) given by 

■?i, 7 (?) = Va,^{q)L^ C {q) + E E jO( N a \y)v a >,yLyj(jq Fa ,^)§ a ,^.^(jq Fa ,^,q) . (A21) 

a', 7' j=±l 



Inserting Eq. (|A14| ) in Eq. (|A21|) we obtain Eq. ©. 
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APPENDIX B: STATIC MASSES FOR THE HEAVY PSEUDOPARTICLES 



The static mass m* is defined in Ref. [|T3[ as 



2t drj a ^{r) / dr 



m 



(2vrp ai7 (r)) 2 



2ti] an (r) (2w dp a ^ (r ) / dr) 



(27rp Qi7 (r))3 



(Bl) 



where the functions 2tr] an {r) and 2np an (r) are defined in Ref. |T5[ and r° is W® (qFa,-/) 



which represents Q, r Cj7 , and r S)7 , 

After some straightforward algebra, the general expressions lead to the following 
simple expressions for 1/m* 

1 _4te 2 /(l + « 2 7 2 ) 3/2 + A2, 7 



and 



1 A 77 -, - M - A 77 



7 >0 



(B2) 



s,7 ■' 1 -s,7+2 



In Eqs. pi) and 



m s,7 (^£,7+1 — As i7 — A^ j7+2 ) 2 

the functions and A£, j:r read 

9f q ,o dq v Qfi (q)R { al(q) 



7 > 



(B3) 



A!!, 



^ [i + (<!>(<?)/*) 2 ] 



and 



A p 



(B4) 



a , 7tx 1 + 



(B5) 



with 



(B6) 



In the limit of fully polarized ferromagnetism, these expressions lead to the following 
closed-form expressions for the static masses 

I tn I i\ + 2[?7i )7 ] M7sin(2n7r) 



m c, 7 %Vy] 2 V sfl + [m 7 ] 2 M 2 + sin 2 (mr)^ 



(B7) 



tTT 



1 m(7 + 1) sin(2ri7r) 

\ vV / l + [n(7 + l)] 2 2 b?2,7+i] [u(7 + l)] 2 + sin 2 (n7r)^ 



^BS) 
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where 



771,3. = tan~ 1 [cot(ri7r) 7 ==== ] , (B9) 



and r)2, x = it/2 — 771^. We remark that the static masses of the 0,7 pseudoparticles are, in 
general, negative. The static masses of the a, pseudoparticles have been studied in Ref. 
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APPENDIX C: CHARGE AND SPIN STIFFNESSES AT ZERO TEMPERATURE 



In this Appendix we show that the direct use of Eq. p7|) leads to the stiffness expressions 
(135) - (137) of Ref. §. 

The calculation of the charge and spin stiffnesses P?| ) requires the expansion of Eq. ( p5[ ) 
and of Eqs. (|A7|), (|A"8|), and ( |A9| ) up to second order in cf>. As in the case of the charge and 
spin current, both he charge and spin stiffnesses can be computed from Eq. (0), and we 
obtain one or the other depending on the coupling constants we choose in Eqs. ([A 71) , ( |A8| ), 
and (|A9|) . Expanding the ground-state energy up to second order in 0, we obtain 



d 2 (E /N a 



d(cf>/N a ] 



2 



— f 1FC dq \2tK°'^(q) sm(K(°\q)) + 2t[K^(q)} 2 cos(if (0) (<?))! , (CI) 
2n J-q Fc 1 J 



<f>=0 



where the function K 0, ^(q) is the second derivative of the rapidity function defined by Eq. 
( |A"7| ) in order to (j)/N a at = 0. The functions K 0, ^(q) and K°^(q) can be written as 

K^{q) = d -^^L%{q) , (C2) 
dq 



and 



respectively. The use of Eqs. (|C2|) and ( |C3|) in Eq. ( |Cl|) leads then to 
d 2 (E /N a )_ 

4>=o 



d(<p/N a ) 



l_ 2tsm(Q)[L° c: t(jq Fc )] 2 

27r.^ 27rp c , (g) ' 1 } 
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where the function L^(jq Fc ) is defined in Appendix [A]. The function L^'o^q) obeys the 
following integral equation 

£°'^(g) = ^2 ./ I^.n U</f.s. iiJ 



2tt ^ 27rp s ,o(r S)0 )(l + [sin - jr s , 



2^ 



nJ- gFsfi q dq> l + [sm(K^(q))/u-R%(q')}^ 1 ' 

which was obtained by performing the type of expansions developed in Appendix [A[ More- 
over, L°'^t(g) is given by 

T 0M ( \ _ 1 V- j[cos{Q)L°J{jq Fc )} 2 



2^ x/., , r/n(0). 



4vr 27rp Si0 (r s , )(l + [(/^(g) - jr fl , )/2] 

- - r afi dq' dR{ ^ {q,) L °^ {q,) 



T2> 



2* /-,,„„ dtf l + [(<i(g)-45(g'))/2]2 
JL p W) cos(^(o)(gQ)Lat(g2 

Introducing Eqs. (|C5| ) and ( |Ctj| ) in Eq. (|C|) we obtain, after some algebra, the following 
expression for the (charge and spin) stiffness 

4*1* = J2 v c>0 [L°^(jq Fc , )} 2 + £ ^o[^' C (j^ s ,o)] 2 , (C7) 
j=±i j=±i 

where the functions L 0, ^^(jqF a ,o) axe defined by Eq. ( |A14j ). After some simple algebra, 
expression ( |U7|) can be shown to be the same as expressions (135) - (137) of Ref. 0. 



APPENDIX D: THE ZERO MAGNETIC-FIELD CASE 

For the case of zero magnetic field it is possible to cast the equations for the energy 
bands, phase shifts, and rapidities in a simpler form. After some algebra, the e SjJ (q) band 
(with 7 = 0, 1, 2, ... , oo) and the e Ci7 (g) band with (7 = 1, 2, ... , 00) can at zero magnetic 
field be rewritten as 



2t 



cos(uR ( °l(q)) 

du ' Ti (u) 

00 cosn(u;) 



(Dl) 
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and 



e° cn (q) = ReAty/l - u*(I&%) + * 7 ) 2 - 

roc e ~iu 

- At / du cos(a; J R( }(g))T 1 (a;) 



where Ti(cj) obeys the integral equation 

/oo 
da/Ti(u/)r(u/,u;) 
-oo 

Here the free term and the kernel read 

1 rQ 

TV (a;) = — / <iA;sin(£;) sinfu; sin(£;)/-u) 
2ir J-o 



and 



r / / \ sin((o; - u')r c , Q ) 



(D2) 



(D3) 



(D4) 



(D5) 



respectively. The kernel ( |D5| ) was already obtained in Ref. (see Eq. (A5) of that 
reference). Equation (Dl), together with the fact that in the limit of zero magnetic field the 
width of the s, 7 > momentum pseudo-Brillouin zone vanishes [see Eq. (|2])], shows that 
the s, 7 bands collapse for 7 > and all values of U and n to the point zero. 

In this limit it is also possible to cast the integral equations for the phase shifts, whose 
expressions are given in Ref. | 15| , in the following alternative form 



r c ,o 



$ c ,o;c,o(r, r 1 ) = -B(r - r') + / dr"<S> c>0;Ci0 (r", r')A(r - r") 

-r c ,o 



(D6) 



^v/fcr') = --tan"^^) + T'" dr"® c ^{r", r') A(r - r") 

71 7 ' 



r c ,o 



(D7) 



$ c ,o; S ,7'( r ' r ') = -<Vy' — tan 1 [sinh(7r/2(r - r'))] + 

J-rofi 7 + 1 



(D8) 



7T 7 J-r c , 7T7 1 + (^j 2 



(D9) 
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$ ,(r r') = -6 ,(r - r') - P° ^ WKV) 



$ ,(r r') = - r^ WlV] mil) 



7 



$ s , 7;c ,o(r,r') =0, (D12) 



$ Sj7;c , y (r, r') = ^ f'° ^ ^Tq^'^ > ( D13 ) 
4 J-r c , cosh(7r/2(r" — rj) 



1 fco dr" $ Ci0;Si y(r",r') 



4 y_ rci0 y + 1 cosh(7r/2(r - r")) ' v ; 



and 



i ^^ 5in M,- : -)] e _, 7 , +7 2 _ + e _ 2 „ + _ 

7r Jo a>(l + e^J 



The functions A(r) and 73 (r) are defined as 



A(r) = lr^2?M (D16) 

Wo 1 + e 2|a;| 



and 



B M = T^™> < D17 > 

7r Jo cull + e^ w i; 

respectively. 

At n = 1 we have that Tx(x) = T^x) = Ji(x/u), where Ji(x/u) is the Bessel function 



of order one, and the bands ( pi|) and ( |D2j ) are obtained in closed form. 
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FIGURES 

FIG. 1. The ratio 

m ci/ m * i as function of U at electronic density n — 0.7 and for values of the 
magnetic field h = H/H c = 0.1 (full line), h = 0.5 (dashed line) and h = 0.9 (dashed-dotted line). 
For other electronic densities, the plots follow the same trends as for n = 0.7. 

FIG. 2. The ratio m p c l /m* c l as function of the electronic density n and for values of the magnetic 
field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The onsite Coulomb interaction is (a) U = 1, 
(b) U = 5, and (c) U = 20. 

FIG. 3. The ratio m p cl /m* cl as function of the the magnetic field h and for values of the 
electronic density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite Coulomb interaction 
is (a) U=l,{h)U = 5, and (c) U = 20. 

FIG. 4. The ratio m p cl /m* c l as function of the the magnetic field h and for values of the onsite 
Coulomb interaction U = 1, U = 2, U = 3, C/ = 5, U = 10, and £7 = 20. The electronic density is 

(a) n = 0.5 and (b) n = 0.9. 

FIG. 5. The ratio rn p /m* as function of U, for electronic density n = 0.7, and for values of 
the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. For other electronic densities, 
the plots follow the same trends as for n = 0.7. 

FIG. 6. The ratio m p /m* as function of the electronic density n and for values of the magnetic 
field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The onsite Coulomb interaction is (a) U = 1, 

(b) U = 5, and (c) [/ = 20. 

FIG. 7. The ratio rn p /m^ as function of the magnetic field h and for values of the electronic 
density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite Coulomb interaction is (a) 
U=l,(b)U = 5, and (c) U = 20. 
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FIG. 8. The ratio rn p c Q /m* c Q as function of the the magnetic field h and for values of the onsite 
Coulomb interaction U = 1, U = 2, U = 3, U = 5, U = 10, and U = 20. The electronic density is 
(a) n = 0.5 and (b) U = 0.9. 

FIG. 9. The ratio rn a s \/m* s l as function of U, for electronic density n = 0.7, and for values of 
the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. For other electronic densities, 
the plots follow the same trends as for n = 0.7. 

FIG. 10. The ratio m a s \/m* s l as function of the electronic density n and for values of the 
magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The onsite Coulomb interaction is 
(a) U = l,{h)U = 5, and (c) U = 20. 

FIG. 11. The ratio rn^/m*^ as function of the magnetic field h and for values of the electronic 
density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite Coulomb interaction is (a) 
U = 1, (b) U = 5, and (c) U = 20. 

FIG. 12. The ratio m a s \/m* i as function of the the magnetic field h and for values of the onsite 
Coulomb interaction U = 1, U = 2, U = 3, U = 5, U = 10, and U = 20. The electronic density is 
(a) n = 0.3, (b) n = 0.5, and (c) U = 0.9. 
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TABLES 



H -> H r H -> n 



! 1 

1- 2 



m c, 7 / m c, 7 27(27-^-1) 27(27+5^.^0) 
m e,o/ m c,0 2tsin(7rn T ) 00 
'"«.7/'"a.7 2( 7 +l)(2 7 +2-^ 7 ) 2( 7 +l)(2 7 +2+2ff >7;3i0 ) 2( 7 +l)(2 7 +2+2gi 7;s >0 =g 

Mass ratios in several limits of physical interest. The function ^ 7 is defined as 7/ 7 = 
2/(7r) tan~ 1 [(sin(n7r))/('u[7 + 1])]. The equations for the static masses m* are given in 
Appendix In the case H — > 0, simple expressions for the parameters Ca 7 -a'o> ^q. (ff4|). 
can be obtained from the results of Appendix [D| The ratios m%*/m* and m£ 7 /m* 7 are 
infinite. The dependence on Z7 and n of the parameter £ has been studied in Refs. |5|J5S|. 
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